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Abstract: 

The collider experiments at the Tevatron and the LHC will allow for detailed investiga- 
c^ \ tions of the properties of the top quark. This requires precise predictions of the hadronic 

production of tt pairs and of their subsequent decays. In this Letter we present for the 
reactions pp, pp ^ tt + X ^ i^i'^ +X the first calculation of the dilepton angular distri- 
bution at next-to-leading order (NLO) in the QCD coupling, keeping the full dependence 
on the spins of the intermediate tt state. The angular distribution reflects the degree of 
correlation of the t and f spins which we determine for different choices of t and f spin 
bases. In the case of the Tevatron, the QCD corrections are sizeable, and the distribution 
is quite sensitive to the parton content of the proton. 



PACS numbers: 12.38.Bx, 13.88.-i-e, 14.65.Ha 



* supported by BMBF, contract 05 HTl PAA 4 
^supported by a Heisenberg fellowship of D.F.G. 



The top quark is by far the heaviest fundamental fermion discovered [|l]] to date. It 
is an excellent probe of the fundamental interactions in the high energy regime that will 
be explored by the upgraded Fermilab Tevatron collider and by the CERN large hadron 
collider LHC. It is expected that very large numbers of top quarks will be produced with 
these colliders: eventually about 10^ top quark-antiquark (tt) pairs per year at the Tevatron 
and more than about 10^ ti pairs per year at the LHC. This will make feasible precise 
investigations of the interactions of top quarks. 

Because of their extremely short lifetime top quarks find no time to form hadronic 
bound states: they are highly instable particles whose interactions are governed by short- 
distance dynamics \^. As a consequence the properties of the top quark and antiquark, in 
particular phenomena associated with their spins, are reflected directly in the distributions 
and the corresponding angular correlations of the jets, W bosons, or leptons into which 
the t and f decay. These distributions are determined by the t and f polarizations and 
spin correlations induced by the production mechanism(s). Furthermore they depend on 
the interactions responsible for the top (anti-)quark decay. Hence the analysis of these 
distributions will be an important tool, once large data samples will be available, to obtain 
detailed information about top-quark production and decay. 

For hadronic pair production the spin correlations of ti pairs were studied to leading 
order in the coupling a^ of Quantum Chromodynamics (QCD) in ref. [0, ^. In particular 
it was analyzed which spin bases are most suitable for the investigation of tt spin corre- 
lations induced by the strong interactions^. There exists also an extensive literature, for 
example [^ and references therein, on how to exploit top-quark spin phenomena at hadron 
colliders in the search for new interactions. The work which we report in this Letter serves 
the purpose of putting predictions of tt spin correlations within the standard model of par- 
ticle physics (SM) on firmer grounds. We analyze the hadronic production of tt pairs 
and their subsequent decays, keeping the full information on the spin configuration of the 
tt state. We extend the existing results by taking into account the next-to-leading order 
(NLO) QCD corrections in the production and the decay of the ff pairs. More specifically 
we consider the channels where both t and t decay semileptonically, 

pp,pp^tt + X^i+£'-+X, (1) 

{i = e,iJ,x), and we predict the following double leptonic distribution at NLO in the 
coupling tti-: 

I d^o 1 , 

= -(l-FBicos0++B2COS0_-Ccos0+cos0_) , (2) 

O(icos0+acos0_ 4 

with o being the cross section for the channel under consideration. In Eq. (0) 0+ (0 ) 
denotes the angle between the direction of flight of the lepton £+ (i'^) in the t (t) rest 
frameQ and a reference direction a (b). The directions a, b can be chosen arbitrarily. 

'a first attempt to measure tt spin correlations with a very small data sample was made in ref. p]. 

^ We define the rest frame of the t{t) quark by a rotation-free Lorentz boost from the center-of-mass 
(cm.) frame of the initial partons that produce the tt pair. If one defines the t{t) rest frame by a boost from 
the hadronic cm. frame, it will differ from our choice by a Wigner rotation. 



Different choices will yield different values for the coefficients Bi,2 and C. The physical 
interpretation of these coefficients is well known [^ ^]: The coefficient C in Eq. (^ 
reflects spin correlations of the tt intermediate state. A more detailed discussion will 
be given below Eq. @. For our choices of the directions a and b (cf. Eq. (§)) QCD 
interactions yield vanishing coefficients Bi, B20. 

In principle one could measure the angular distribution of every possible decay prod- 
uct of the top (anti-)quark. In the SM, where the main top-quark decay modes are 
t -^ bW -^ bqq',bi\£, the most powerful analyzers of the polarization of the top quark 
are the charged leptons, or the jets that originate from quarks of weak isospin —1/2 pro- 
duced by the decay of the W boson. Here we restrict ourselves to the double leptonic 
distribution. 

To predict the "dilepton -i- jets" distribution (Q) at NLO accuracy we have to consider 
the following parton subprocesses: 

gg,qqJUbbi+i'-VeVe,, (3) 

gg, qq ^ bbi+i'-\i\e + g, (4) 

g + q{q)^ bbi+i'-Veyt +q{q). (5) 

At the Tevatron the cross section is dominated by quark-antiquark annihilation while at 
the LHC gluon-gluon fusion is predicted to be the dominant production process. 

In view of the fact that the total width Yt of the top quark is much smaller than its mass 
nit (^t/i^t = 0(1%)), one may analyze the above reactions using the so-called leading 
pole approximation [|^]. This amounts to expanding the amplitudes of Eqs. (|3|) - (^ around 
the poles of the unstable t and f quarks. Only the leading term of this expansion, i.e., 
the residue of the double poles is kept here. The radiative corrections to the respective 
lowest-order amplitudes can be classified into so-called factorizable and non-factorizable 
corrections. We take into account the factorizable corrections to the above reactions for 
which the squared matrix element 5Vf is of the form |fAf p »« Tr[p7?p]. Here R denotes the 
respective spin density matrix for the production of on-shell tt pairs, and p (p) is the t (f) 
decay density matrix. 

To obtain a theoretical prediction for the distribution in Eq. @ at NLO accuracy 
we use our recent results [0] on the tt production spin-density matrices at NLO QCD. 
These results extend previous calculations [||] of the differential tt cross section with 
spins summed over and allow the calculation of the cross section for a specific spin con- 
figuration of the tt state. In particular, the quantization axes can be chosen arbitrarily. 

The decay density matrix p (p) required for computing @ describes the normalized 
angular distribution of the decay of a polarized t{t) quark into i^{i^) + anything in the 
rest frame of the t{i) quark. The matrix p has the form 2pot/„ = ( 11 + K_|_ a • q+)a'a where 
q+ describes the direction of flight of i^ in the rest frame of the t quark and a, denote the 
Pauli matrices. The decay matrix p is obtained from p by replacing q+ by — q and K+ 

^ This is due to the parity invariance of QCD. 



by K_. The factor k+ (k_) signifies the top-spin analyzing power of the charged lepton. 
It is equal to one to lowest order in the SM, that is, for V —A charged currents. Its value 
including the order a^ corrections can be extracted from the results of [ |T0| ] and turns out 
to be very close to one: K+ = K_ = 1 — O-OlStti. Using the general expressions for p, 
p and the fact that the factorizable contributions are of the form Tr[pi?p] one obtains the 
following formula for the correlation coefficient C in Eq. @: 

C = 4K+K_((a-s,)(b-sr)), (6) 

where Sf,Sf denote the t and t spin operators. The expectation value in Eq. @ is defined 
with respect to the matrix elements for the hadronic production of tiX. It is related to the 
more familiar double spin asymmetries 

AI(^ . Vh . ^\ - N(TT)+N(ii)-N(n)-N(iT) ,_. 

^^ ^ ^^^ ^ ^^^"N(TT)+Nai)+N(n)+NaT)' ^ ^ 

where A'^(tT) ^tc. denote the number of Z^? pairs with t and t spin parallel - or anti-parallel 
- to a and b, respectively. From Eq. (^ one can see that the axes a, b introduced through 
the angles Q± in Eq. (g) can be interpreted as quantization axes of the intermediate ti state 
within our approximation. Eq. @ generalizes the lowest-order results of ^, 0] and holds 
for factorizable contributions to all orders in the QCD coupling^. 
For definiteness we consider here the following spin bases: 

a = kf, b = kf (helicity basis), 

a = p, b = p, (beam basis), (8) 

a = d/, b = df (off-diagonal basis). 

Here k^ (k^) denotes the direction of flight of the t{i) quark in the parton c.m.s., and p is 
the unit vector along one of the hadronic beams in the laboratory frame. Furthermore df 
is the axis constructed in ref . [Q] with respect to which the spins of t and f produced by qq 
annihilation are 100 % correlated^ to leading order in a^. (For gg -^ tt one can show that 
no spin basis with this property exists.) 

Table [l] contains our results^ for C at leading and next-to-leading order in a^ using 
the parton distribution functions (PDF) CTEQ5L (LO) and CTEQ5M (NLO) of [0. 
(These numbers and the results given below were obtained by integrating over the full 
phase phase. Results with cuts included will be given elsewhere [jOJ].) For pp collisions 
at a/* = 2 TeV the helicity basis is not the best choice because the t, t quarks are only 
moderately relativistic in this case. Table [l] shows that the dilepton spin correlations at 
the Tevatron are large both in the off-diagonal and in the beam basis. In fact they are 

^The non-factorizable NLO QCD corrections were calculated for gg and qq initial states in ref. [pl|. We 
expect with these results that the effect of these corrections on the dileptonic angular correlations is small. 

^ We use the definitions for df and df given in ref. [R|. In particular, dj = dr at LO. The sign of CoS. at 
LO is therefore opposite to that of [B| . 

^We use the MS factorization scheme, a^ is defined to be the five-flavour MS coupling, and m, is defined 
in the on-shell scheme. 



almost identical. The QCD corrections decrease the LO results for these correlations by 
about 10%. Since the gg initial state dominates ;^? production with pp collisions at y/s = 
14 TeV the beam and off-diagonal bases are no longer useful. Here the helicity basis is a 
good choice and gives a spin correlation of about 30%. In this case the QCD corrections 
are small. The large difference between the LO and NLO results for the correlation in 
the beam basis at the LHC is due to an almost complete cancellation of the contributions 
from the qq and gg initial state at LO. 





pp aty^^ 


= 2 TeV 


pp at y/s = 


= 14 TeV 


LO 


NLO 


LO 


NLO 


Chel. 


-0.456 


-0.389 


0.305 


0.311 


Cbeam 


0.910 


0.806 


-0.005 


-0.072 


Coff. 


0.918 


0.813 


-0.027 


-0.089 



Table 1: Coefficient C ofEq. (^ to leading (LO) and next-to-leading order (NLO) in a^ 
for the spin bases ofEq. (^. Theparton distribution functions of^T^ were used choosing 
the renormalization scale jjr equal to the factorization scale jjf = mi = 175 GeV. 
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Figure 1 : Dependence of Of Cbeam ^^ LO (dashed line) and at NLO (solid line) on /u 
I^R = Pf for pp collisions at y/s = 2 TeV, with PDF of j^^. 



We now discuss the uncertainties of our predictions. It is well known that the inclusion 
of the QCD corrections reduces the dependence of the tt cross section Ot on the renor- 
malization and factorization scales significantly. The same is true for the product o^C. 



"400 



200 



10 



dashed line: LO 
full line: NLO 



10 



Figure 2: Dependence o/OfChei. at LO (dashed line) and at NLO (solid line) on jj ■ 
jup for pp collisions at ^ = 14 TeV, with PDF of [^D^. 
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In Figs. |T]and ^we demonstrate this with OfCbeam and o,Chei. evaluated at Tevatron and 
LHC energies, respectively, as functions of /j/mt, where ij = /jr= /jf. The corresponding 
figure for OfCoff. is almost identical to Fig. |l]. 

To leading order in a^ the coefficient C depends only on the factorization scale /Uf, 
while at NLO it depends on both scales /j^ and /jf. Table ^ shows our NLO results for 
the three choices /jr= /jf = mt/2,mt,2mt, again using the PDF of [[T^. An extension of 
this work, which is however beyond the scope of this Letter, would be the resummation 
of Sudakov-type logarithms at the next-to-leading logarithmic level. This was performed 



in ref . Qlq ] for the total cross section Ot and it stabilizes the predictions for O/ with respect 
to variations of /jr and fjf. 



IJR=IJF 


ppat^^l TeV 


ppat^= 14 TeV 


*-hel. *-beam *-off. 


Chel. 


mt/2 


-0.364 0.774 0.779 


0.278 


nit 


-0.389 0.806 0.813 


0.311 


Init 


-0.407 0.829 0.836 


0.331 



Table 2: Dependence of the correlation coefficients, computed with the PDF of j^^, on 
f^ = l^R = f^F at NLO. 



In Table |] we compare results for C using different sets of PDF. In the case of pp col- 
lisions at ^/^ = 2 TeV, the spread of the results is larger than the scale uncertainty given 
in Table |^. To a considerable extent this is due to an interesting feature of C, namely 
the qq and gg initial states contribute to C with opposite signs. Therefore the spin corre- 
lations are quite sensitive to the relative weights of qq and gg initiated tt events. These 
weights depend in particular on the chosen set of PDF. For example, one finds the fol- 
lowing individual NLO contributions for the helicity, beam, and off-diagonal correlation 
at the upgraded Tevatron: for the GRV98 (MRST98) PDF C^fi = -0.443 (-0.486), 



Cik = +0.124 (+0.075), Cll^ = +0.802 (+0.879), C^^ = -0.068 (-0.042), and 



C^l = +0.810 (+0.889), C^ff = -0.073 (-0.044). This suggests that accurate measure- 
ments of the dilepton distribution @, using different spin bases, at the upgraded Tevatron 
may provide additional constraints in the continuing effort to improve the knowledge of 
the PDF 



PDF 


pp at y/s = 2 TeV 


pp at ^= 14 TeV 


^hel. ^beam ^off. 


Chel. 


GRV98 


-0.325 0.734 0.739 


0.332 


CTEQ5 


-0.389 0.806 0.813 


0.311 


MRST98 


-0.417 0.838 0.846 


0.315 



Table 3: Correlation coefficients Chei., Cbeam. <^nd Coff. at NLO for i^r = i^f = nit (^nd 
different sets of parton distribution functions: GRV98 /|7^, CTEQ5 /[I?|/, and MRST98 
(c-g) 0. 



Finally we have studied the dependence of the C coefficients on the top quark mass. 
For this we have used again the CTEQ5 PDF and set fj = nit. In the case of pp collisions 
at ^/s = 2 TeV, a variation of mt from 170 to 180 GeV changes Chei. from —0.378 to 
-0.397, Cbeam from 0.790 to 0.817, and Coff. from 0.797 to 0.822. At LHC energies, 
Chel. changes by less than a percent. 

The extension of our results to the "lepton-hjets" and "all jets" decay channels [ [r3| ] is 
straightforward. The "lepton-hjets" channels should be particularly useful for detecting 
tt spin correlations: although one looses top-spin analyzing power one gains in statistics 
and the experimental reconstruction of the t and t rest frames may also be facilitated. 

In conclusion we have analyzed, at next-to-leading order in a^, the hadronic produc- 
tion of tt quarks in a general spin configuration and have computed the dileptonic angular 
correlation coefficients C that reflect the degree of correlation between the t and t spins. 
Our results for the Tevatron show that the scale and in particular the PDF uncertainties in 
the prediction of the dileptonic angular distribution must be reduced before tt spin cor- 
relations can be used in a meaningful way to search for relatively small effects of new 
interactions that are, for example, not distinguished by violating parity or CP invariance. 



Our results may also be useful to learn more about the parton distributions in the proton 
at high energies. For pp collisions at ^/s =14 TeV the theoretical uncertainties in the 
prediction of this distribution are smaller and one may adopt the optimistic view that at 
the time the LHC will be turned on further theoretical progress will have turned top quark 
spin correlations into a precision tool for the analysis of tt events. 
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